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NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


MEMORANDUM 1-17“59A 

ESTIMATION OF STATIC LONGITUDINAL STABILITY OF AIRCRAFT 
CONFIGURATIONS AT HIGH MACH NUMBERS AND AT ANGLES 
OF ATTACK BETWEEN 0° AND ±l80° 

By Duane W. Dugan 

SUMMARY 


The possibility of obtaining useful estimates of the static longitudinal 
stability of aircraft flying at high supersonic Mach numbers at angles of 
attack between 0° and ±l80° is explored. Existing theories , empirical formu- 
las, and graphical procedures are employed to estimate the normal-force and 
pitching-moment characteristics of an example airplane configuration consist- 
ing of an ogive -cylinder body, trapezoidal wing, and cruciform trapezoidal 
tail. Existing wind-tunnel data for this configuration at a Mach number of 
6.86 provide an evaluation of the estimates up to an angle of attack of 35° ■ 
Evaluation at higher angles of attack is afforded by data obtained from wind- 
tunnel tests made with the same configuration at angles of attack between 
30° and 150° at five Mach numbers between 2.5 and 3*55* Over the ranges of 
Mach numbers and angles of attack investigated, predictions of normal force 
and center-of-pressure locations for the configuration considered agree well 
with those obtained experimentally, particularly at the higher Mach numbers. 


INTRODUCTION 


The design of an aircraft which is required to leave and re-enter the 
atmosphere poses a number of problems. For example, it is conceivable that 
re-entry into the atmosphere may be made advertently or inadvertently at 
any angle of attack between 0° and ±l8o°. An angle of attack of the order 
of 30° or bo° may be desirable in order to decelerate rapidly without undue 
heating, or to develop sufficient lift to skip out of the atmosphere. 

Angles of attack between 0° and ±l8o° may also be encountered by air-to-air 
missiles launched by high-speed aircraft. 
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The purpose of this report is to explore the possibility of obtaining 
useful estimates of the static longitudinal stability of aircraft configu- 
rations at high supersonic Mach numbers and at angles of attack between 0° 
and ±l80°. To do this, existing theories, empirical formulas , and graphical 
procedures are employed to predict the normal-force and pitching-moment 
characteristics of a given airplane configuration over a range of angles of 
attack from 0° to ±l80°. The configuration for this study is shown in fig- 
ure 1. Wind-tunnel data for this airplane up to an angle of attack of 35° 
at a Mach number of 6.86 are available for limited evaluation of the 
estimates. To obtain an indication of the accuracy of the predictions at 
higher angles of attack, a model of the same configuration is tested at 
angles of attack from 30° to 150° at each of five Mach numbers between 2.5 
and 3-55- 

It should be noted that the specific theories and procedures employed 
in this report are somewhat arbitrarily selected, and other choices are not 
precluded. 


ANALYSIS 


The aerodynamic characteristics of an aircraft configuration may be 
calculated by summing the aerodynamic characteristics of the isolated 
individual components and the effects due to the interference of one part 
upon the aerodynamic characteristics of the others. Thus the equations 

% = % + c % + C Kt + AC % £ W j + AC Kb ( T j + ACj^ (B J + AC %(b ) + AC % ^ } + AC^^ 

( 1 ) 

Cm = C mB + C mw +C mT +AC mB ^ w j +ACn^^ T ^ +AC mw ^ B j +AC mT ^ B j +AC n]ir ^ w j +AC mw ^ T j 

( 2 ) 


provide a convenient framework for obtaining the static longitudinal 
stability of aircraft configurations consisting of various combinations 
of body, wing, and tail. The subscripts B, W, and T pertain to the 
body, exposed wing, and exposed tail, respectively; the terms AC^g/y ) , 


ACjj w ( b )> etc., represent the interference effects of the exposed wing 

upon the body, of the body upon the exposed wing, etc. The notation used 
herein is given in appendix A. 


The next several sections are devoted to developing procedures for 
estimating over a range of angles of attack from 0° to ±l8o° and for high 
Mach numbers the following: (l) normal forces and pitching moments of a 

fuselage consisting of a pointed or rounded nose section followed by a 
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cylindrical portion; (2) normal forces and pitching moments of an isolated 
wing (or tail); and ( 3 ) effects of interference of one part of a 'body-wing- 
tail configuration upon the normal forces and pitching moments of the 
others . 

It should be clear that modification or extension of the procedures 
used here might be required in assessing the aerodynamic characteristics 
of configurations having features dissimilar to those of the one considered. 
Obviously, a number of the steps and calculations included in the present 
study would not be required in the case of a more simple conf iguration 
such as an all -wing conf iguration or a pointed body of revolution with 
stabilizing fins. 

Forces and Moments on the Body 

At high Mach numbers, the flow about a body of revolution can be 
approximated by the Newtonian concept which assumes that the individual 
particles comprising the fluid medium possess only translational velocities 
in the direction of flow and lose all their momentum in a direction normal 
to those surfaces on which they impinge. Based on this impact theory, the 
coefficient of pressure at any point on surfaces exposed directly to the 
flow is given by (see ref. l) 

Cp = 2(sin 0 cos a- sin a cos G sin p) 2 ( 3 ) 

where, as shown in sketch (a) below, 

Q angle made by surface of body with body axis 
a angle of attack of body axis 

P polar angle of any point on body surface, measured from positive xy 

plane and positive for counterclockwise directions when viewed from 
rear 



Sketch (a) 
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The value of the pressure coefficient on shielded surfaces Cp^ 

■would he zero in the impact theory. Consideration of gas dynamics indicates 
that the value of C p will vary between zero and -2 /rH» depending upon 

the free-stream Mach number, shape of the body, and angle of attack. ^Accord 
ing to reference 2, the lowest pressure in the flow about bodies is, from, 
experimental evidence at subsonic and low supersonic Mach numbers, approxi- 
mately 0.3 of the free-stream static pressure. In wind-tunnel tests (re . 3J 
of a circular cylinder at angles of attack from 0 to 90 and at a Mach. num- 
ber of 6.86 the lowest pressure measured on the lee surfaces of the cylinder 
was approximately 0.3^ of the free-stream static pressure. For a value of 
7 equal to 1.4, then, the largest negative value 1 of Cp u appears to be of 
the order of -l/M^ 2 . The mean value of the pressure coefficient over the 
shielded surfaces will thus be very small for large Mach numbers and can be 
considered as zero without significant effects on the calculated values of 
normal force and pitching moments . 

For any portion of the fuselage shown in sketch (a), the normal force 
is given by 


N = - 2 <3oo 


dx 1 


'Pu V 

C-nr sin pdp + / C p r sin pdp J 
J Pu 


'Pu 




^length ''• u - «/ 2 

or, with C p given by equation (3), and with Cp^ assumed zero, 

dC w 2r \(_ *Y ~ . 1 

—21 = — sin 2a cos 2 0 
dx S 


P u + IJtan 9 + ± cos p u ( cot a tan 2 0 + 2 tan a 


OO 


where S is a reference area. Likewise, the pitching moment, taken about 
the centroid of area of the base of the nose section for convenience, is 
given by 


M = 


'length! 


Ijj - x ) - r tan 0 




cL xf f U C p r sin pdp + £ Cp u r sin pdp 


or 


■'■Values of C p negatively greater than -l/Mo, 2 have been observed 

at supersonic speeds in tests of wings at moderate. and large angles of 
attack; however, the magnitude of the coefficient is generally closer to 
-l/M^ 2 than to -2/ /M^ 2 • . 
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dC m 2r 


sin 2a cos ^0 { - x j - r tan 0 ^ p u + 2 A 311 0 + 


cos p u (cot a tan 2 0 + 2 tan a' 


where l is a reference length. 

Equations ( 4 ) and ( 5 ) may he integrated analytically or numerically 
to give coefficients of normal force and pitching moment for any arbi- 
trarily shaped body of revolution. Calculations for a conical nose are 
particularly simple, since 0 = constant = 0 V , the semivertex angle of 
the cone, and r is a linear function of x. For positive angles of 
attack less than 0 V , p u = it/2 and all surfaces are exposed to the flow, 
(i.e., Cp is defined for all values of p). For 0 V < cl < Jt - 0 V , P or ~ 
tions of the cone are shielded from the flow, and p u = sin" 1 (tan 0 v /tan a). 
For it - 0 V < a < jt, no surfaces of the cone are exposed to the flow and 
p u = -at/2 . In the above ranges of angles of attack the corresponding 
normal-force and pitching-moment coefficients based upon the area of the 
base and upon the length of the cone are 


r sin 2a ; 


0 < a < Q v 


cone it 


— cos 2 0 v sin 2a p u + — + 


^ cos P u (cot a tan 0 V + 2 cot 0 v tan a)J ; 


9 V < cl < k - 6^ 


Ct\J — 0 \ Jt “ 0 ~\t ^ CL ^ J( 

^ cone v — — 


— (l - 2 tan 2 0y )cos 2 0 v sin 2a ; 0 < a < 0 A 


(l -2 tan 2 0 v )cos 2 0 v sin 2a Pu + ^ + 


— cos P u (cot a tan 0 V + 2 cot 0 v tan a) ; 0 V < a < it - 0 


0 J jt — 0-y ^ a < C jt 
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From the above results, the location of the center of pressure of 
the conical nose is independent of angle of attack; its distance from the 
vertex is given by 

= \ (l + tan 2 0 v ) (8) 

l N 8 

Calculations for other nose shapes are more involved since r, 9, 
and p u are all functions of the lengthwise variable x. Inasmuch as 
the body-wing- tail configuration under consideration (fig. 1) has an 
ogival nose section, the integrations indicated in equations (4) and (5) 
are carried out for the circular-arc ogive with the following results: 


Cw =| C 2 sin 2cd (l - if 2 ) 2 - 2*Jl - q 2 + \ (l + sec a) - 

^ogive i 2 


Jr- csc 2 a sec aJl - ti 2 
3* 1 


(1 + k 2 )E(k) 



0 < a < 6 V 

(9a) 


C N • = TT C 2 sin 2aj(l - r] 2 )[§ + sin -1 7 — 

"ogive 3 it (_ X 2 tan ajl - t} 2 


2rtJl - T ] 2 + y 


3(l + rj 2 )sin 2 a - t] 2 + b 


. - 1 ^ it 8 

3 sec a sin x H — + 3 o “ T 
J sin a 2 3 


esc a sec a cos cp + 

(1 +k 2 )E(cp,k) 


csc 2 a sec cxa/ 1 - r) 2 


k ,2 F(cp,k) 


fl v < a < j - 0, 


(9b) 


Cw . = % C 2 sin 2a-| J (l+sec a) - 2*JT- t] 2 

"ogive 3 I 2 


_8 

3 « 


csc 2 a sec awl-rj 2 


(1+ k 2 )E(k) - k ,2 K(k) 


> Tt “ Q-y ^ CC ^ Ti 

(9c) 


and 
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C 2 / 

Cm * = Tz~r si - T] 2 sin 2a 


Jm ogive 15 t) 

5t]( 3 + 2r| 2 ) 


40 o So / 0 . o \ 0 

— a csc^a - -r— (3 - Sinhalese 2a - 
Jt 3* 


s/l - rj2 - s j_ n 2a + 15 ^sin~ 1 rj - + 


(3 + sin 2 a)csc 2a - ^ cos a cot 2 a cosh~ 1 sec a + — 1 Fi L : 
* v jt 56* x f ' 


0 < a < 


(10a) 


Cm . = 7 : — \/l - ti 2 sin 2a( 3 24 csc 2 a[a - tan™ 1 (tan a cos cp ) ] - 

^ogive 9^3 v 1 


16(3 - sin 2 a)csc 2a - 6t] (3 + t] 2 )sin~ 1 
2(3 + T] 2 ) (l - 3 sin 2 a) 


9 + 7 sin^a + 


tan a*s/l - t] 2 
esc a sec a cos cp - 3^(3 + 2 t] 2 ) *Vl - t] 2 - 


9 cos a cot 2 a(^cosh~ 1 sec a-cosh“ 1 sec a\/l- r] 2 ^ + 

rt 3 \ 

6(3 + sin 2 a)csc 2a + 9 2 s ^- n_ lr l + 280 ^ 2 ) 9 ^v < a £ n “ 


(10b) 


q 2 j 

Cm . = ™ — s/l - *n 2 sin 2a-s6(3 + sin 2 a)csc 2a - 9 r sin 2a + — (jt - a) - 

^give 9^1 1 4 2 


9 cos a cot 2 a cosh -1 sec a - 8 


2 ( 3 - sin 2 a)csc 2a + 


3 csc 2 a(jt - a) 


3 


n/ 1 ” + 2 Qq ^ < a < 7t 


(ioc; 


In the above equations 

f N Z N 
11 C “ R 


R 

d 
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R radius of generating arc of ogive 

fjy fineness ratio of nose — 


k sin a 


9 


sin 


T) 

sin a 




4w 4 m 

qnd 2 ’ qjtd 2 lu 


and F-l, F 2 , and F 3 are given in appendix B. 


If f N = C = l/2 , we have the hemispherical nose, and the above 
expressions for and C m become 


and 


C N 


hemisphere 


1 

4 


(sin 2a + 2 sin a) ; 


0 < a < at 


Cm, • t_ =0; 0 < a < jt 

m hemi sphere — — 

In passing, it is interesting to compare the initial normal-force 
curve slopes of the cone and ogive for a given fineness ratio fjj. From 
the above results for the ogive, 



4c 2 

^ (cos^0y *" 4 COS Qy + 3 ) 


2 j" 1 + 32f^ 4 

3 L(l + 4f w 2 ) 2 _ 


and from equation (6a) for the cone. 



- 2 cos 2 0 v 


8f lf 

(l+4f N 2 ) 


The ratio of these slopes is 

( c NcJ 

( c NcJ 


cone 


ogive 


12f N 2 (l + 4f N 2 ) 
l +8f n 2 + 48f N 4 
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which is unity for fjj = l /2 and oo, is greater than unity for l /2 < % < oo, 

and has a maximum value of only 1.06 when fjj = I/ 2 J 1 + n/ 17 = 0.777- For 
angles of attack near zero, therefore, the normal force of an ogive can he 
adequately obtained from the simple equation for that of the cone of equal 
fineness ratio. For fineness ratios of unity or larger, calculations show 
that the differences in Cjj for the cone and ogive are negligible at 
angles of attack up to 9 V for the cone. At angles of attack somewhat 
less than 9 V of the cone, the curves of Cjy versus a for the two nose 
shapes cross, so that at a = 9 V , Cn for the ogive exceeds that for the 
cone . 


Wo explicit expression for the location of the center of pressure can 
be given for the ogive as in the case of the cone; computations have shown, 
however, that for small angles of attack the center of pressure of the 
ogive is nearer the vertex than that of the cone of equal fineness ratio 
and moves rearward with increasing angle of attack. 

The contributions of the cylindrical portion of the fuselage to the 
normal force and pitching moment of the complete body when calculated from 
equations (A) and (5) are 


and 


C N, 


cyl 


16 L 

3* d 


s in 2 a ; 


0 < a < 3 t 




cyl 


8 L L . 2 
•x — -y t — sm^a ; 

3* i In 


0 < a < jt 


(id 


( 12 ) 


where L is the length of the cylinder, and the coefficients are based 
on the same reference area and length as in the case of the nose section 
above . 

Equations (ll) and (12) do not include the effects of centrifugal 
forces in the flow about the cylinder, nor the effects of viscous forces 
(skin friction). These effects tend to cancel each other as shown by 
the results of tests on cylinders reported in reference 3 . Equation (ll) 
predicted normal forces in excellent agreement with the experimental 
force data at small to moderate angles of attack and overestimated such 
forces by only 5 percent near 90°. On the other hand, inclusion of the 
effects of centrifugal forces in the impact theory (ref. 1 ) gave values 
lower than those obtained experimentally by approximately 10 percent at 
the lower angles of attack and 5 percent near 90 °. 
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Forces and Moments on Isolated Wing and Tail 


Forces and moments on the isolated ving and horizontal tail are 
estimated over the range of angles of attack from 0° to l8o° hy a combina- 
tion of two-dimensional shock-expansion theory (corrected for tip effects) 
and the Newtonian impact theory modified as discussed below. 

For angles of attack not too near nor beyond the angle at which the 
bow wave of the wing or tail is detached , shock-expansion theory provides 
a good approximation to the lifting pressures in the absence of tip 
effects. The coefficient of normal force based on the exposed area of 
the wing is first calculated from two-dimensional shock-expansion theory 
for angles of attack up to that for shock-wave detachment. Next, if the 
reduced aspect ratio (3A is so small that tip effects become appreciable, 
the two-dimensional results are multiplied by the ratio of the lift-curve 
slope of the finite wing to that of a wing of infinite span, both values 
being obtained from linearized supersonic theory. Values of lift-curve 
slopes for finite wings may be obtained readily from such sources as 
reference 4 . The values of thus obtained are plotted as a function 

of the angle of attack. Forces on a tail in the wake of a wing and vice 
versa are estimated subsequently in the section entitled "interference 
Effects." 

To the writer* s knowledge, no methods exist for successfully predict- 
ing the lifting pressures on finite wings at angles of attack beyond that 
for shock detachment. The Newtonian impact theory, so useful in estimating 
normal forces on bodies of revolution for all angles of attack and even at 
comparatively low supersonic Mach numbers, is very inadequate in the case 
of planar surfaces. Even at very large Mach numbers the impact theory 
grossly underestimates the lifting pressures on an infinite flat plate at 
small to moderate angles of attack; on the other hand, the theory over- 
estimates the force on a plate placed normal to the free stream according 
to experimental evidence. In order to arrive at some logical basis for 
estimating the normal forces on a finite wing at angles of attack beyond 
which the shock-expansion method is inapplicable, the following approach 
is adopted. 

First, the assumption is made that the variation of normal force on 
a flat surface at angles of attack in the immediate neighborhood of 90° 
is that given by the Newtonian concept; that is, the force varies as the 
second power of the sine of the angle of attack. Next, the impact theory 
is modified by replacing the factor 2 in the equation C]\j = 2 sin 2 a by 
a factor in accord with experimental results for flat surfaces placed 
normal to the free supersonic stream. The third step is to fair an 
interpolating curve from the curve obtained above from considerations of 
shock-expansion theory (with any corrections necessary for tip effects) 
to the curve based on the modified Newtonian impact theory. This latter 
procedure will be discussed more in detail subsequently when an application 
to the wing of the example body -wing-tail configuration is made. 
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The choice of a factor to replace the factor 2 in the impact theory 
is made from the following considerations. Experimental measurements 
(refs. 5 > 6 , and 7 > and unpublished data) of the pressure over the face 
of right circular cylinders with axes alined with the free supersonic 
stream indicate that the pressures remain nearly constant at the value 
given by gasdynamics for stagnation pressure p s over approximately half 
the distance from the center of the face to the edges, and thereafter 
decrease to a value p so corresponding to sonic velocity at the edge. 

The variation of the ratio of local pressures p to stagnation pressure 
p s with fractional distance £ from center to edge is shown in figure 2 
for faces of cylinders of various diameters at several Mach numbers. 

Also shown in the figure is a plot of the equation 


V_ 

^s 



Pso 
Ps . 


1(1 -£ 4 ) 


1/3 


which appears to fit the data within a few percent. The corresponding 
expression for the coefficient of pressure on the windward face is 




7K 


f -i 

-^oo 


~ c p. 


SO 


L C P S 


+ 1 - 


^Pso 

: Ps 


(1-i 4 ) 


1/3 


where 

I 

C Pso 

C Ps 

C Pso 

cp 

Pso 


ratio of distance from geometric center to distance from center to 
edge 

coefficient of pressure corresponding to sonic velocity 
coefficient of pressure corresponding to stagnation conditions 

-i 


so 


hin- — 


0.528 for 7 = 7/5 in air 


The coefficient of normal force for the windward face of the cylinder 
C jj based on the area of the face is then 




CpjS 
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Cii, 



0.842 + 0.158 



A comparison of the last equation with integrated experimental pressure 
distributions from several sources is presented in figure 3 . Inasmuch as 
the distribution of pressure over a flat surface normal to the free super- 
sonic stream is apparently determined primarily by the acceleration of the 
air from rest at the geometric center to sonic velocity at an edge, the 
equations given above for the flat face of circular plan form may be 
expected to serve as reasonable estimates for a variety of plan forms . 

Now consider the pressures on the lee surfaces of flat plate shapes 
normal to the free stream. If viscosity is ignored, a Prandtl-Meyer 
expansion from sonic velocity would, of course, result in zero pressure 
over the lee surfaces, giving a pressure coefficient = -2//M . As 

noted in an earlier section, there exists a limit pressure coefficient 
in air which has been found experimentally to have a value close to 
-l/Moo 2 . The latter value is adopted here as the average pressure coeffi- 
cient on the lee surfaces of flat shapes normal to supersonic flow. The 
total maximum normal-force coefficient is then 


C N, 


'max 



0.842 + 0.158 




(13) 


Figure 4 shows a comparison of the above equation with experimental data 
obtained in ballistic firings and in wind-tunnel tests over a wide range 
of Mach numbers. These data are taken from references 8 , 9, 40, and 11, 
and from unpublished results of tests in the Ames 10- by l4-inch wind 
tunnel. The data obtained in the Ames 10- by l4-inch wind-tunnel tests 
at nominal Mach numbers 3, 4, 5, an( i 6 indicate that the plan form of a 
flat wing at an angle of attack of 90 ° has but little effect on the total 
streamwise force. Pressure measurements made on the lee surfaces of the 
wings in the same tests substantiate the choice of -l/M » 2 as the average 
value of Cp. From figure 4 it is concluded that the empirical equation 
derived here satisfactorily predicts normal-force coefficients for wings 
of various plan forms at a = 90° not only at large Mach numbers but at 
Mach numbers as low as 2. No reasons can be suggested here for the 
occurrence of a minimum value for C N between Mach numbers 3 and 4 as 
indicated by the data; in any case, the discrepancies between predicted 
and experimental values of C K are not so large, percentagewise, as to 
preclude the use of the empirical equation. The variation with Mach 
number of the stagnation pressure coefficient used in the equation is 
shown in figure 4 for comparison. 

For the present purpose, the variation of the normal-force coeffi- 
cient of an isolated wing or tail at angles of attack in the immediate 
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neighborhood, of 90° will be estimated as C I -j rnay sin 2 a, with given 

by equation ( 13 ). There remains the problem of estimating the normal 
forces in the angle -of -attack range between angles of attack near shock 
detachment and those near 90°. As noted earlier, graphical interpolation 
is used in the absence of theory or experimental data. Experimentally 
obtained variations of normal force with angle of attack up to and beyond 
shock detachment may serve as guides. For example, figure 5 presents the 
normal -force characteristics of a rectangular wing of aspect ratio 3 a t 
angles of attack well beyond shock detachment for three different Mach 
nu mb ers (ref. 12). Also shown in the figure are the modified impact 
values at a = 90° for each of the Mach numbers, together with interpolat- 
ing curves. It is noted that at angles of attack well below those for 
shock detachment, the normal-force curves become increasingly concave with 
increasing Mach number as would be expected from predictions of two- 
dimensional shock-expansion theory; at angles of attack in the neighbor- 
hood of shock detachment, however, the curves are nearly linear, and 
become definitely convex at still larger angles. The above observations 
may also be made in the case of the normal-force characteristics of a 
triangular wing of aspect ratio 4 presented for two Mach numbers in fig- 
ure 6 (ref. 13). The modified impact-theory values and interpolating 
curves are again drawn. Also shown in the figure is the variation of 
normal-force coefficient with angle of attack at Moo = 6.86 for the wing 
of the example configuration as calculated from two-dimensional shock- 
expansion theory. Although the profile of this wing has a rounded leading 
edge and a blunt trailing edge, for the purpose of calculation it was 
assumed that the profile was a symmetrical wedge -slab -wedge with sharp 
leading and trailing edges. The angle of attack at which shock detachment 
theoretically occurs is indicated in each case. The above angles were 
calculated according to the method presented in reference l4, which takes 
into account both the thickness of the profile and the leading-edge sweep 
of the wing. With the curves of the delta wing serving as qualitative 
guides, an interpolating curve for the example wing was drawn as shown . 
Values of normal-force coefficient for the isolated wing at angles of 
attack between 0 ° and l 80 ° can then be obtained from the composite curve. 

In the case of the tail of the example configuration, the wedge pro- 
file, although rounded, was assumed to be sharp-edged for purposes of 
calculation by means of two-dimensional shock-expansion theory. For rear- 
ward flight, because of the high Mach number, the shock wave was assumed 
to be effectively attached to the lower surface until the angle of flow 
deflection was equal to the angle calculated from two-dimensional shock- 
wave theory for shock detachment. The modified impact theory was based 
upon the angle made by the lower surface of the wedge with the free 
stream rather than the angle of attack. 

To estimate the contributions of the isolated wing or horizontal 
tail to the pitching moment of the configuration, the eenter-of-pressure 
location is calculated by shock-expansion theory for angles of attack at 



lb 


which the bow wave remains attached, 2 and thereafter extrapolated to the 
value given by impact theory at a = 90°. In figure 7 the variation with 
angle of attack of the distance of the center of pressure from the leading 
edge of the mean aerodynamic chord of the exposed wing is plotted up to 
the angle at which the bow wave detaches. As in calculating normal forces 
above, the wing profile is assumed to be sharp-edged and symmetrical with 
respect to the midchord. Therefore, the center of pressure of the wing 
at a = 90 ° is at the centroid of area of the exposed wing, or at midchord 
of the mean aerodynamic chord. Since the tail of the example configura- 
tion has a wedge profile, the center of pressure is assumed to be at the 
centroid of area of the exposed tail at all angles of attack. 


Interference Effects 


I | In this section estimates are 

c f made for the effects upon normal 

I force and pitching moments of (l) 

I % the wing and tail upon the body; 

/ — -= ^ (2) of the body upon the wing and 

x s /3y yl I ' tail; and ( 3 ) of the wing upon the 

h } I tail (and vice versa) for angles of 

I I attack between 0° and l8o°. In 

/ / most cases only first-order effects 

/ / h ~ Cr are considered; that is, the fact 

I h* x* £y + c f that the influence of the wing 

I upon the body, for example, would 

/ / affect the net influence of the 

j 1 body upon the wing which in turn 

/ / y would affect the influence of the 

H.- wing upon the body, and so on, is 

^ d ^ neglected here since the extra 

^ d ^ labor involved would not be justi- 

M fied in the present method of esti- 

i mation (cf. ref. 15)* In assessing 

Sketch (b) the influence of the horizontal 

tail upon the body in the case of 
the body -wing-tail configuration, however, the normal forces on the tail 
in the wake of the wing rather than in the free stream are used since no 

extra labor is involved and the effect upon the pitching moment of the 

configuration may be appreciable in some applications. 


Body in presence of wing or tail .- In order to estimate the inter- 
ference effects of the wing and tail upon the body, linearized supersonic 
theory is employed. The method of approximation outlined in reference 1 6 
is adapted to the present problem. Sketch (b) above shows the essential 

2 These calculations can be used subsequently in determining the 
characteristics of the wake of the wing. 
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feature of the procedure in the case of a trapezoidal wing or tail of 
large span mounted on a cylindrical body. Note that for simplicity the 
cylindrical body is approximated by a plane surface of the same plan form; 
thus, the Mach cones which define the region of influence of the wing 
intersect the body in straight lines. The ratio of the normal force on 
the body due to the wing to that of the isolated exposed wing is given by 
(for 0° < a < 90°) 
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where 


P n/Hx , 2 "I 

m cotangent of sweep angle of leading edge of wing 


h length of cylindrical portion of body from juncture of body and wing 
leading edge to base or pd + c r , whichever is smaller 


c r root chord of exposed wing 

c t 

A e taper ratio of exposed wing, 37 
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s semispan of wing-body combination 

normal-force curve slope of exposed isolated wing from linearized 
supersonic theory (see ref. k ) 

Equation (15) is in its most general form. If d 2 - h/p (the case shown 
in sketch (b)), the equation is somewhat simplified since the last two 
terms vanish. Likewise, if h = pd + c r , then d x = d 2 = d, and the 
equation reduces to that given in reference 1 6. 

The point at which this additional body normal force may be considered 
to act, x/c r measured from the juncture of the leading edge of the wing 
and body, is found from 
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where Mg(w)> "the moment of the induced normal force on the body about the 
juncture of the leading edge of the wing and the body, is calculated to be 
generally 
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and is obtained as above from linearized supersonic theory, 

cx v 
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If, as shown in sketch (b), d .2 = h/p, the last four terms in equa- 
tion (17) vanish. If h = pd + c r , then d x = d 2 = d, and equation ( 17 ) 
reduces to that given in reference 1 6 . 

The interference effects of the wing (or tail) upon the body in 
reversed flight (90° < a < l 80 °) are found in the same manner as above with 
the proper values of m and h. 


The interference normal forces and moments will be given by 
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, coefficient of normal force of the exposed wing as 


estimated above from shock-expansion theory and modified impact 
theory 


Cjfn 


Nrp 




, coefficient of normal force of the isolated tail in the 


case of wingless configurations as estimated above from shock- 
expansion theory and modified impact theory 


L %(w) 

S 

x cg 
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Nt ill } coefficient of normal force for the tail in the wake of 

1;he wing for 0° < a < 90 ° in the case of body-wing-tail 
configurations as estimated below 

reference area 

distance of moment reference from nose of fuselage 

distance from nose of fuselage of the cp of the additional 
normal force on the body due to the wing 
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distance from nose of fuselage of the cp of the additional 
nor ma] force on the body due to the horizontal tail 

reference length 

Obviously, the above method cannot be used for angles of attack for 
which the shock wave is detached from the wing or tail. At an single of 
attack of 90°, impact theory implies no interference of the wing or tail 
upon the body; at finite high Mach numbers for which the detached bow wave 
at angles of attack near 90 clings closely to the wing or tail, the assump- 
tion that the above interference is negligibly small appears reasonable from 
consideration of gasdynamics. In order to arrive at some estimate of the 
additional forces and moments on the body due to the wing and tail over the 
ranges of angles of attack from 0° to l 80 °, resort to graphical interpolation 
is again made as shown in figures 8 through 13 > where the procedure is 
applied to the example configuration at a Mach number of 6.86. In these 
figures, the curves for the aerodynamic characteristics of the body includ- 
ing wing or tail interference are faired smoothly into the corresponding 
curves of the body alone. The normal-force and pitching-moment coefficients 
of the body alone were calculated by adding the contributions of the ogival 
nose and of the cylindrical portion of the body obtained as discussed in the 
preceding section. The interference effects were obtained as indicated by 
equations ( 15 ) through (21), with |3C Na = ^ as for a two-dimensional super- 
sonic wing. If it is assumed that the above interference is given closely 
by linearized supersonic theory for angles of attack not too near that for 
detachment of the bow wave and decreases gradually and monotonically toward 
zero as the angle of attack approaches 90 °, then the difference between any 
two interpolating curves which might be employed to join the pairs of curves 
representing the aerodynamic characteristics of the body with and without 
interference is not likely to cause more than slight differences in final 
results for the complete configuration. 

Wing and tail in presence of body .- As pointed out in reference 17, 
at very high supersonic speeds the normal force on a cylinder-wing combina- 
tion (neglecting wing interference on the body) is approximately equal to 
the normal force of the wing alone if that portion of the wing covered by 
the cylinder is included. Here the covered portion of the wing is con- 
sidered to be an extension of the root chord of the exposed wing through 
the cylinder in a diametral plane. The additional normal force on the 
exposed wing is, therefore, the difference between the normal force which 
would be carried by the portion of the wing covered by the body and that 
carried by the cylinder having a length equal to the root chord of the 
exposed wing. 

From consideration of impact theory alone (Moq-^oo) , the normal force 
on a flat plate wing having an area Cjyd is 

% c = 2q O0 sin 2 ac rw d 


X AB(T) 
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where 

Cjy root chord of exposed wing 

d diameter of cylindrical body 

Likewise, the normal force of the cylindrical body of length Cp^ and 
diameter d is 


N cyl = | ( 3co sin2cw: i V d 

The interference normal force on the exposed wing due to the presence of 
the body is then 


^W(B) = 3 ^ sin ^iy d 
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( 22 ) 
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%(B) 2 L w " C yl 

, coefficient of normal force of cylindrical portion of body 


cy! <3 c^B 

with length L as estimated in a preceding section 
In the same way, the interference normal force on the tail is estimated as 
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According to reference 17, these additional forces are concentrated 
in narrow regions adjacent to the wing-body and tail-body junctures. The 
pitching moments resulting from the interference normal forces on wing 
and tail due to the presence of the body can, therefore, be approximated 
by assuming that these forces act at the midchord position of the respec- 
tive root chords of the exposed wing and tail. 

The above estimates of interference effects are considered applicable 
in the present case over the entire range of angles of attack from 0° to l8o°. 

Tail in wake of wing.- In order to estimate the loads on a horizontal 
tail in the wake of the wing of a configui rtion flying at high supersonic 
speeds, the angle of flow, Mach number, and dynamic pressure of the wake 
must be determined. The procedure adopted here to obtain the above charac- 
teristics of the wake is described below as applied to the body-wing-tail 
configuration of figure 1. 
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Sketch (c) shows the salient features of inviscid supersonic flow past 
a two-dimensional flat-plate wing of chord length AB at an angle of attack 
a. In the region ABC the fluid is compressed and turned to flow parallel 
to the lower surface; in the region CBF it next undergoes expansion and a 
turning back toward free-stream direction, and crosses the final expansion 
ray BF at an angle u) with respect to the chord plane. In general, the 
angle w is larger than the angle of attack a. Beyond the point of inter- 
section with the ray BF, the fluid eventually returns to free-stream direc- 
tion. Above the wing, in the region ABG, the fluid is expanded and turned 
to flow parallel to the upper surface; as it crosses the terminal shock 
wave BHG, the fluid undergoes recompression and is turned toward the free- 
stream direction again, although its angle with respect to the wing chord 
as it crosses BHG is generally larger than a. Right at the trailing 
edge B, the fluid adjacent to the upper and lower surfaces of the wing 
leaves at the angle u) and follows the streamline BE. This streamline 
is often termed a "slip line" since there are discontinuities in velocity, 
Mach number, temperature, etc., (but not in pressure) across it. Reflec- 
tions such as JK of the expansion waves between AH and AG from the 
shock wave BHG serve to return the streamline from the trailing edge to 
the free-stream direction at some point E as shown. 


Although upwash occurs right at the trailing edge of the wing, the 
flow in the region GBC aft of the trailing edge may be characterized as 
either upwash or downwash depending upon the location in the region, upon 
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the free-steam Mach number, and upon the angle of attack of the wing. The 
other characteristics of the flow such as Mach number and dynamic pressure 
also depend upon the same parameters. The dynamic pressures of the wake 
decrease from comparatively high values in regions adjacent to BC to 
relatively low values above BF, while the Mach number of the wake increases 
with increasing distance from BC toward BE. Hence, as has been pointed 
out elsewhere (e.g., ref. l8), the effectiveness of a horizontal tail 
placed in that region of the wake above the chord plane of the wing will 
be generally less than that of one located below. In what follows, the 
horizontal tail is assumed to be located on or below the chord plane of 
the wing. 

The pertinent characteristics of the flow behind the trailing edge of 
a two-dimensional wing at angles of attack less than that for shock detach- 
ment are estimated here on the basis of shock-expansion theory. For a 
given free -stream Mach number, the characteristics of the flow along any 
ray of the expansion fan CBF (sketch (d)) may be found conveniently by a 


A 



graphical procedure as follows : For a given angle of attack ex, plot a 
curve of 0 versus Av with the latter as the independent variable from 
the relation 


where 


0 = Av + p 1 - b 2 
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Av v 2 -v x 

v 1 ,v 2 angle through which air at M = 1 must turn in Prandtl-Meyer flow 
to accelerate to M^Mg 

u ,u Mach angles, sin -1 , sin ' 1 77 - 

1 2 Mg 

The value of Av for the flow across any ray is then read from the curve 
at the value of 0 for the ray in question. Since M x , Vx/v^f and v x are 
known"* for each angle of attack a, the value of v 2 is determined and 
hence M 2 and P 2 /Pi may be found from Prandtl-Meyer relationships (most 
conveniently from tables of compressible flow such as are included in 
ref. 19)- The ratio of the dynamic pressure of the flow at any point along 
a ray to that of the free-stream may be computed from 

h. _ dkf AA AA 

%> VW \pj \vj 

For a tail whose chord plane is alined with one of the rays in the expansion 
fan CBF, the direction, Mach number, and dynamic pressure will be constant 
along the chord. For other orientations of the tail, the above characteri- 
istics of the flow will vary along the chord of the tail; in such cases cal- 
culations for the forces at any spanwise station of the tail may be simpli- 
fied somewhat without appreciable error by using the average of the flow 
characteristics at the leading and trailing edges, or by assuming that the 
flow calculated at the midchord is constant over the entire chord of the 
tail. For tails with zero dihedral angle, the spanwise distribution of 
pressure will be constant (excluding tip effects); for tails with positive 
or negative dihedral, the pressures will diminish or increase toward the 
tip, respectively, and graphical or numerical integration may be necessary 
to determine the total normal force on the tail. 

At some Mach numbers and at small angles of attack, the tail may be 
outside the expansion fan CBF. In this case the flow approaching the tail 
will have completed its expansion and turning, and is then subject to 
recompression and a deflection in the opposite direction (toward free 
stream) through the influence of reflected disturbances from the flow over 
the upper surface of the wing such as the wave KL shown in sketch (c). 
However, for small angles of attack the above effects of the flow will not 
be large and in the few cases for which this situation exists the flow may 
be regarded as having the characteristics it had at the terminal expansion 
ray BF. The latter ray may be determined graphically from the stipulation 
that the common direction of the flow at the trailing edge coming from both 
lower and upper surfaces must be such that the recompression of the air 
passing over the upper surface and the expansion (in general) of the air 
from the lower surface result in the same value for the static pressure. 

^These three quantities can be readily obtained from previous 
calculations for the center of pressure of the wing. 
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Certain aspects of the flow behind a wing at high supersonic speeds 
have been discussed in references l8 and 20. As noted in such references, 
the nature of the flow in the wake of a wing is determined to a consider- 
able extent by the thickness of the wing; consequently, the direction, 
static pressure p x , and Mach number M 1 used in the above procedure are 
calculated for the actual profile of the wing rather than for a flat plate. 

The characteristics of the flow behind the wing of the configuration 
of figure 1 in the region of the tail were calculated according to the 
above procedure and are presented in figure l4. The profile of the wing 
was considered to be a wedge -slab -wedge for these calculations. The odd 
variation with angle of attack noted for the dynamic pressure in the wake 
is due to the combination of nonlinear relationships among the various 
parameters involved in shock-expansion theory. The final static pressure 
of a high-speed gas which is first subjected to a nonisentropic compression 
and then expanded isentropically depends in a nonlinear fashion upon both 
the extent of the compression and the amount of expansion. At angles of 
attack near that for shock detachment (40° in the present case) the increas- 
ing losses in total head across the oblique shock wave from the leading 
edge of the wing combined with increasing expansions in the wake were 
apparently sufficient to produce the downward trend in the variation of 
q T noted in figure l4 at angles of attack near 4o°. 

Since the calculations for the normal force on the tail in the wake 
of the wing cannot extend beyond that angle of attack for which the shock 
wave detaches from the leading edge of either the wing or the tail, some 
means of estimating the force on the tail for the larger angles must be 
devised. At angles of attack near 90°, the wing will not appreciably 
affect the tail behind it. This fact suggests the use of graphical inter- 
polation to approximate the force on the tail at large angles of attack 
as illustrated in figure 15. As indicated, the variation with angle of 
attack of C N for the tail in the wake of the wing determined from shock- 
expansion theory is plotted up to the angle of attack for bow wave detach- 
ment. Also plotted in the figure is the same variation for the isolated 
tail in the free stream computed as explained in an earlier section from 
a combination of shock-expansion theory and modified impact theory. A 
third curve is then drawn to effect a transition from one computed curve 
to the other. The composite curve then represents the sum of the normal 
force of the tail and the interference normal force due to the wing and 
will be related to the term Cu T + ACn t ( w ) of equation (l) in final 
calculations . 

The contributions of the tail to the pitching moment of the body- 
wing-tail configuration will be based on the normal force of the tail in 
the wake of the wing and on the assumption that the center of pressure of 
the tail remains at the centroid of area of the exposed tail. 
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Wing in vake of tail . - The effects of the tail upon the aerodynamic 
characteristics of the wing (disregarding effects of the tail tips) m ay 
be estimated in much the same way as above. Inasmuch as the span of the 
tail is generally considerably less than that of the wing, the effect of 
the wake may extend over only a portion of the wing; consequently, the 
net effect on the normal force of the configuration may be comparatively 
small and the effect upon pitching moment is likely to be negligible for 
body-wing-tail configuration in which the wing is located close to the 
center of gravity. 

The effect of the flow around the tail tips is confined within narrow 
Mach cones originating at the tail tips, so that at large Mach numbers the 
portion of the wing affected by such flow is likely to be comparatively 
small. In addition, the portion of the wing affected by the tail tip 
vortices will be divided into nearly equal regions of upwash and downwash, 
so that the net effect upon normal force and pitching moment of the complete 
configuration will generally be negligibly small. 

Since for the present body-wing-tail configuration the effects of the 
wing upon the tail were not large (fig. 13 ), the effects of the tail upon 
the wing in reversed flight were judged to be too small to justify the 
calculations involved and were thus ignored. 


COMPARISON BETWEEN PREDICTED AND EXPERIMENTAL 
RESULTS FOR THE AIRPLANE CONFIGURATION 


An evaluation of the foregoing estimates of normal -force and pitching- 
moment characteristics of the airplane configuration selected for this study 
is now made. Unfortunately, experimental data available at any given Mach 
number for comparison with the estimates do not extend over the complete 
range of angle of attack from 0° to l80°. The airplane configuration con- 
sidered here has been tested in the NASA Langley 11-inch blow-down wind 
tunnel at a Mach number of 6.86 up to angles of attack of 35° • (Results 
up to a = 28° are published in ref. 21.) Figures l6, 17, and l8 give the 
comparison between predicted and experimental variations of normal force, 
pitching moment, and center-of -pressure location with angle of attack at 
Moo = 6.86, not only for the complete airplane but also for the body alone 
and the body -wing and body-tail combinations. As far as the experimental 
data permit comparison, the estimated results are in satisfactory agreement 
with the wind-tunnel tests in each case. A comparison between the results 
both with and without the inclusion of interference effects (fig. 19 ) 
indicates that inclusion of such effects should be made at the lower hyper- 
sonic Mach numbers such as 6.86. 

To obtain an indication of the accuracy of the estimates at higher 
angles of attack the following procedure was followed. Available wind- 
tunnel facilities enabled testing of a model of the example configuration 
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over a range of nominal Mach numbers between 2.5 and 3*5- Accordingly, a 
model was constructed and tested at five Mach numbers in the given range at 
angles of attack between 30° and 150° in the test section of the Ames 8- by 
7-foot Unitary Plan wind tunnel. The tests are described in appendix C 
of this report. The basic data obtained in the tests, normal-force coef 
ficients and center-of -pressure locations, are presented in figures 20 
and 21. The procedures described in the preceding sections were then 
applied to the complete configuration to give the estimated variation of 
normal— force coefficient and center-of —pressure location over a range of 
Mach numbers from 2.5 to 100. Predicted and experimental results (from 
faired curves of figs. 20 and 21) for a number of angles of attack between 
30° and 150° are shown as a function of l/k^ in figures 22 and 23 . In 
general, the agreement between the estimated and experimental variation 
of normal force and center-of -pressure location with Mach number is good, 
not only qualitatively but quantitatively. However, at angles of attack 
near 90 ° ( 75 ° to 105 °), the experimental variation of normal force differs 
somewhat from that predicted. Whereas a monotonic decrease of normal force 
coefficient with Mach number is predicted at all angles of attack, the 
present data suggest that the coefficient may have a minimum value Q in the 
neighborhood of M^ = 3*5 for angles of attack between 75 and 105 • It 
may be recalled that, similarly, an apparent minimum in the normal-force 
coefficient of wings at a = 90° was indicated between Mach numbers 3 an d 

4 by the data presented in figure 4. As a matter of interest, averaged, 
experimental values of Cjj from figure 4 were used for the wing and tail 
to calculate the normal-force coefficient of the complete example configure 
tion at an angle of attack of 90° • The results are shown as the dotted 
curve in figure 22. Although the agreement between the calculated and 
experimental variation of normal force is generally improved by using 
experimental values for the wing and tail in this case, it is obvious 

that other factors are involved. Quantitatively, the discrepancies between 
the estimated and the present experimental values of for the configura 

tion at angles of attack between 75° and 105° are generally less than 

5 percent except at angles near 100 ° where they are of the order of 10 per- 
cent. There is reason to believe, moreover, that the accuracy of the 
estimates of normal force for this configuration at angles of attack near 
90° may be expected to improve with increasing Mach number beyond a Mach 
number of approximately 3 * 5- 

Another phenomenon observed in the present tests which is not fully 
understood may be seen in figure 20. The abrupt decrease in normal force 
as the angle of attack increases from 90° to 100° is certainly not pre- 
dicted as shown in figure 24 where the predicted variation of normal-force 
coefficient with angle of attack at a Mach number of 2.5 is compared with 
experimental data. The corresponding discrepancies between predicted and 
experimental pitching-moment curves are shown in figure 25 . That the 
nature of the flow about the model changes as the attitude of the model 
changes from forward flight to rearward flight can be seen in a typical 

of schlieren pictures taken during the present tests and presented 
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in figure 26. It is also possible that the shrouded sting and balance 
case attached to the lee side of the model in the tests may have contrib- 
uted to disturbances in the flow over portions of the model in the subject 
angle-of -attack range more than anticipated. Whatever the causes, the 
data in figure 18 suggest that the distortion noted in the normal-force 
curves diminishes with increasing Mach number. 


CONCLUDING REMARKS 


Existing theories, empirical formulas, and graphical procedures have 
been employed in an attempt to predict normal-force and pitching-moment 
characteristics of a given airplane configuration at high supersonic Mach 
numbers and at angles of attack ranging from 0° to l8o°. Availability of 
wind-tunnel data at a Mach number of 6.86 provided an evaluation of the 
predictions up to an angle of attack of 35 ° • To obtain an indication of 
the accuracy of the estimates at higher angles of attack, wind-tunnel 
tests of the same configuration were made at angles of attack between 30 ° 
and 150° over a range of Mach numbers from 2.5 to 3«55* 

From a comparison between predicted and experimentally obtained values 
of normal -force coefficients and center-of -pressure locations over the 
range of Mach numbers and angles of attack investigated, it is concluded 
that for Mach numbers of the order of 5 and larger, close estimates of 
these characteristics are possible up to at least 150°. At Mach numbers 
between 2.5 and 3 - 55 , the estimates agree well, in general, with the 
results of the wind-tunnel tests at angles of attack from 30 ° to 150 °. 

From the manner in which the estimates are obtained here at large 
angles of attack, application of the present procedures does not appear 
practical for Mach numbers much lower than 2.5, particularly for angles 
of attack larger than 90 °. 

In assessing aerodynamic characteristics of conf igurations having 
features dissimilar to those of the one considered, it should be clear 
that modification or extension of the procedures used in this report might 
be necessary. Obviously, a number of the steps and calculations included 
in the present study would not be required in the case of a more simple 
configuration such as an all -wing configuration or a pointed body of 
revolution with stabilizing fins. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Oct. 17 , 1958 
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APPENDIX A 
NOTATION 


mean aerodynamic chord of entire wing or tail 
root chord of exposed wing or tail 

maximum diameter of nose section; diameter of fuselage 
fineness ratio of body of revolution 
modulus of elliptic integrals (sin a) 
reference length 

cotangent of angle of sweepback of leading edge of wing or tail 

pressure at a given point on surface of body, wing, or tail 

static pressure in free stream 

. . vM 2 p 

dynamic pressure, — - — 

local radius of body of revolution 

semispan of wing or tail of complete configuration 

longitudinal, lateral, and vertical coordinates of a right-hand 
Cartesian coordinate system 

P 

caliber of ogive of revolution, ^ 


coefficient of pitching moment, based on reference area S and 
on reference length l 

coefficient of normal force, based on reference area S 


coefficient of pressure, 


P " P 


<Jco 




stagnation pressure coefficient, 
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E(cp,k) 

E(k) 

F(cp,k) 

K(k) 

L 

K B(W) 

M 

N 

R 

S 

v 

¥ 00 

a 

P 


7 


n 

ACm B(W) 


AC 


“B(T) 


ACm T(B) 


ACm T(W) 


coefficient of pressure corresponding to sonic velocity 

coefficient of pressure in expansion regions of flow 

incomplete elliptic integral of second kind 

complete elliptic integral of second kind 

incomplete elliptic integral of first kind 

complete elliptic integral of first kind 

length of cylindrical portion of fuselage 

ratio of normal force on body due to presence of wing to normal 
force of exposed wing alone 

pitching moment ; Mach number 

normal force 

radius of generating arc of circular- arc ogive 
reference area 
velocity of free stream 
angle of attack of body axis 

radial angle of a point on surface of body of revolution; 
cotangent of Mach angle 

7 

ratio of specific heats, — for air 

liq 

ratio of length to radius of ogive, 

additional pitching-moment coefficient of body due to presence 
of wing 

additional pitching-moment coefficient of body due to presence 
of tail 

additional pitching-moment coefficient of tail due to presence 
of body 

difference between pitching-moment coefficient of tail in wake 
of wing and that of isolated tail 
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ACm W(B) 

AC mw(T) 

A %(W) 

a %(t) 

ac %(b) 

AC %(w) 

ac %(b) 

ac %(t) 

A 

0 

<?v 

9 

B 

c 

eg 

cp 

cyl 


additional pitching-moment coefficient of wing due to presence 
of body 

difference between pitching-moment coefficient of wing in wake 
of tail and that of isolated wing 

additional normal-force coefficient of body due to presence of 
wing 

additional normal-force coefficient of body due to presence of 
tail 

additional normal-force coefficient of tail due to presence of 
body 

difference between normal-force coefficient of tail in wake of 
wing and that of isolated tail 

additional normal-force coefficient of wing due to presence of 
body 

difference between normal-force coefficient of wing in wake of 
tail and that of isolated wing 

taper ratio of wing or tail 

angle made by surface of body of revolution with body axis 


semivertex angle of nose section 
argument of elliptic integral , sin 



Subscripts 


body or fuselage 
cone 

moment reference 
center of pressure 
cylinder 

exposed wing or tail 



lover surface of ving or body 


nose section of body 
horizontal tail 
ving 


free-stream conditions 
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APPENDIX B 


DERIVATIONS OF F x , F 2 , F 3 

The expressions F x , F 2 , and F 3 of equation (10) in the text are the 
result of applying the method of integration in series to one intiactahle 
integral. The integral in question is 


I 



tan a 


_t 

■s/l - 1 2 


dt 



2 


t 4 

2-k 


lutf.liiil... Y sln -x 1= 

2*4*6 2*4*6*8 /\ tan ajl - t 2 


where 0 < t 2 < Z n 2 /R 2 < 1* For most applications, only four terms of the 
series expansion above are necessary for good accuracy. Based upon^the 
above approximate evaluation of the integral, the terms F x > F 2 , and F 3 
are 


F 1 = 2 s i n a 


l 680 - (280 + 42 sin^a+15 sin^a)sin 2 a 


1343 a - sin a cos a 


337+ (38 + 8 sin 2 a)sin 2 a 


Fp = T] 


1680 - (280 +42r) 2 + 15 ^ 4 )ti 2 
\l si^ 2 ^ - 


sin 


- 1 


, , , _n -V sm^a - t 

1343( a - tan 


cos a 


- sm a cos a 


tan <xn /1 - rj 2 

337 + (38 + 8 sin 2 a)sin 2 a 


cos a 


337 + (8 sin 2 a + 38 + 4r] 2 )sin 2 a + (3”n 2 + 


\l sin 2 a - T ] 2 


and 


^ 3t . , 

F 3 = 2 sin °°l 


(280+42 sin 2 a + 15 sin 4 a)sin 2 a - 1680 


1343 (at - a) - sin a cos a 


337+ (38+ 8 sin 2 a)sin 2 a 
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APPENDIX C 


EXPERIMENTAL INVESTIGATION OF THE STATIC LONGITUDINAL 
STABILITY OF THE AIRPLANE CONFIGURATION AT ANGLES 
OF ATTACK BETWEEN 30° AND 150° AND AT MACH 
NUMBERS BETWEEN 2.48 AND 3-55 


As was stated in the text of this report, wind-tunnel data for the 
body -wing-tail configuration considered extend to only approximately 35° 
angle of attack at a Mach number of 6.86. To obtain an indication of the 
accuracy of the estimates at larger angles, the present investigation was 
conducted over a range of angle of attack from 30° to 150°. Inasmuch as 
the Mach number range of the available testing facility does not extend 
beyond 3.55, measurements were made at five different Mach numbers between 
2.48 and 3.55 in order to determine the trend of normal force and center- 
of-pressure travel with Mach number. 


Apparatus 


The tests were conducted in the 8 - by 7 -foot test section of the Ames 
Unitary Plan wind tunnel. This wind tunnel is a continuous -flow, single- 
return type. The stagnation pressure can be held constant at values 
between 5 and 56 inches of mercury, and the Mach number can be varied dur- 
ing operation between nominal values of 2.4 and 3 . 5 . The sting support 
has a range of ±14.65° in the pitch plane. Further details of this facility 
may be found in reference 22 . 

The model was a 4 / 3 -scale replica of the airplane model tested in the 
Langley facilities. Figure 1 and table I give the essential features of 
the Langley model. For testing at angles of attack between 30° and 60°, 
the model was mounted from the downstream end on a 45° bent sting; for 
angles between 75° and 105°, it was supported by a straight sting attached 
dorsally at the midlength of the fuselage; for angles between 120° and 150°, 
the model was mounted as for the previous angle range except that a 45° 
bent sting was used in place of the straight sting. Each model sting was 
supported from a cylindrical balance case housing a six-component strain- 
gage-type balance. The balance itself was supported from the tunnel sting 
support. Both the balance case and model sting were shrouded to prevent 
aerodynamic forces other than those acting on the model from affecting the 
balance readings. A fouling indication system was provided to insure that 
the data would not be affected by mechanical interference between the 
shroud and model, sting or balance case. 
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In order to take into account any increments in the angle of attack 
of the model due to aerodynamic loads during the tests , a cathetometer was 
used to view the model through a window in the test section of the wind 
tunnel. By means of this instrument , the vertical distance between two 
points marked near either extremity of a midmeridian of the model body 
could be determined, and hence the actual angle of attack could be 
calculated. 


Testing Procedure 


For each of the three combinations of model and stings described above, 
readings of the fore-and-aft normal-force and of the axial-force balance 
gages were recorded by a strain-gage printer at each of several angular 
settings of the tunnel sting support for each of the five Mach numbers 
2.48, 2.77, 3-07, 3*30, and 3 . 55. Several readings were repeated. The 
Reynolds number, based upon the mean aerodynamic chord of the model wing, 
varied from 240,000 to 510,000 during the tests. 


Reduction of Data 


The balance forces were resolved into the normal force on the model 
and into the position of the center of pressure. The model normal force 
was then reduced to standard NASA coefficient form based on the area of 
the entire wing of the model. The location of the center of pressure was 
referenced with respect to the nose of the model and made dimensionless 
in terms of the body length. 


Precision of Results 


Because of the manner in which the balance forces were resolved in 
each of the three model- sting-balance configurations, and because of the 
wide range in the magnitude of the forces on the model, the accuracy with 
which the values of C N and x cp /z B could be obtained varied from one 
range of angles of attack to another. On the basis of the sensitivity of 
the instrumentation employed and upon the methods used to resolve the 
balance forces, as well as from consideration of the repeatability of 
results, the maximum errors in the measured quantities in the ranges of 
angles of attack noted are estimated to be 



3 ^ 


C N 


x cp 


a 


( 


30° 

< a < 60 ° 

±0.010 

75° 

< a < 105° 

±0.012 

120° 

k 

< a < 150° 

±0.024 

30° 

< a < 60 ° 

± 0.027 

75° 

< a, < 105° 

±0.012 

120° 

< a < 150° 

± 0.025 

30° 

< a < 60 ° 

±0.15° 

\ 75° 

< a < 105° 

±0.20° 

120° 

< a < 150° 

±0.15° 


v 


It is not known just what effects the sting shroud had on the flow 
over the lee side of the model j however ^ it is Believed such effects were 
negligible with regard to total normal forces and relatively small as far 
as the position of the center of pressure was concerned. 
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TABLE I.- GEOMETRIC CHARACTERISTICS OF LANGLEY MODEL (REF. 21 ) 


Wing 

Area (including area submerged in fuselage ), sq in 6.240 

Area (exposed), sq in 4.402 

Mean aerodynamic chord (entire wing), in 1.713 

Mean aerodynamic chord (exposed wing), in 1.456 

Span, in 4.330 

Root chord (entire wing), in 2.530 

Root chord (exposed wing) , in 2.133 

Tip chord, in 0.35^- 

Aspect ratio (entire wing) 3-000 

Sweep of leading edge, deg 38*830 

Sweep of c/4 line, deg 29*000 

Incidence at fuselage center line, deg 0 

Dihedral, deg 0 

Geometric twist, deg 0 

Horizontal and vertical tails 

Area (including area submerged in fuselage), sq in 2.o6o 

Area (exposed), sq in 1.204 

Span, in 2.690 

Mean aerodynamic chord (entire tail), in 0,853 

Mean aerodynamic chord (exposed tail), in 0.571 

Root chord (entire tail), in 1.2l4 

Root chord (exposed tail), in 0.950 

Tip chord, in 0.317 

Aspect ratio 3*5^0 

Sweep of leading edge, deg 22.63 

Dihedral, deg 0 

Fuselage 

Length , in 7 • 500 

Maximum diameter, in 0.790 

Fineness ratio 9*500 

Base diameter, in 0.790 

Distance from nose to moment reference, in 3*950 

Maximum cross-sectional area, sq in 0.490 

Ogive nose length, in 2.290 

Ogive radius, in 6.850 
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38.83 



Figure 1.- Three-view sketch of example airplane configuration tested in Langley 11-inch 
hypersonic wind tunnel (ref. 21); all dimensions in inches. 









Figure 2.- Pressure distribution over flat faces of circular cylinders 

with axes parallel to free stream. 




kl 



Mach number 

Figure 3 • “ Variation with Mach number of normal force due to pressures 
on windward faces of cylinders with axes parallel to free stream. 



Mach number 


Figure k.- Variation with Mach number of total streamwise force on bodies 
having flat faces of various plan forms; flow normal to flat faces. 







Figure 5 • ~ Normal-force characteristics of a rectangular wing of aspect ratio 3 at angles of 
attack up to and beyond those for shock detachment for three Mach numbers . 




configuration at = 6,86. 
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Figure 7.- Estimation of position of the center of pressure for the wing of the example 

configuration: M„ = 6.86. 
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Figure 8.- Estimation of the normal force on the body in the presence of the wing in the case of 

the example configuration; = 6.86. 









Figure 11.- Estimation of the pitching moment of the body in the presence of the horizontal tail 
in the case of the example body-tail configuration; = 6.86. 




Figure 12.- Estimation of the normal force of the body in the presence of the horizontal tail in 
the case of the example body -wing- tail configuration; M = 6 . 86 . 




the case of the example hody-wing-tail configuration; M = 6.86. 








Angle of attack a, deg 

Figure Flow characteristics in region of tail in the wake of the wing of the example 

configuration; = 6.86. 




Normal-force of coefficient 



Figure 15.- Estimation of the normal force of the tail in the wake of the 
wing in the case of the example configuration; = 6.86. 


Normal-force coef f icient Cn . q 



Angle of attack a, deg 


Figure l6.- Variation of normal-force coefficient with angle of attack for 
the example airplane configuration and its components; = 6.86. 



Normol-force coefficient C 
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Figure 1 6 .- Concluded. 


Pitching- moment coefficient 



Figure 17.- Variation of pitching-moment coefficient with angle of attack 
for the example airplane configuration and its components; = 6.86. 



Distance of 



Angle of attack a ,deg 

Figure l8.- Variation of location of center of pressure with angle of attack 
for example airplane configuration and its components; = 6.86. 



Distance 



Figure l8.- Concluded. 





o 10 20 30 40 50 60 70 80 90 

Angle of attack a ,deg 


Figure 19-“ Comparison "between results obtained by including and neglecting 
interference effects in estimating normal force and pitching moment for 
example airplane configuration; Moo = 6.86. 
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Figure 21.- Experimental variation of center-of -pressure location with angle of attack for example 

airplane configuration at five supersonic Mach numbers. 





configuration with Mach number at a number of angles of attack. 




Figure 22. - Continued. 







(deg) 



Figure 22.- Concluded 






Figure 23-- Estimated and experimental variation of center-of -pressure location with Mach number 
for example airplane configuration at a number of angles of attack. 



Figure 23.- Concluded 
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Angle of attack a, deg 

Figure 24.- Variation of normal-force coefficient of example airplane configuration with angle of 

attack; = 2.5. 






Figure 25.- Variation of pitching -moment coefficient of example airplane configuration with angle 

of attack; = 2.5- 
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a = 105° a = 150° 

Figure 26.- Schlieren photographs of example airplane configuration indicating 
change in nature of flow with attitude of the model; = 2.77. 
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